Nodes of the superconducting gap probed by electronic Raman scattering in 

HgBa2CaCu20Q^s single crystals 
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Pure electronic Raman spectra with no phonon structures superimposed to the electronic con- 
tinuum, are reported, in optimally doped HgBa2CaCu20%+s single crystals (21 = 126 K). As a 
consequence, the spectra in the A\g, Big and -B29 symmetries, including the crucial low energy 
frequency dependence of the electronic scattering, are directly and reliably measured. The -629 and, 
most strikingly, the Big spectra exhibit a strong intrinsic linear term, which suggests that the nodes 
are shifted from the [110] and [110] directions, a result inconsistent with a pure d^2_y2 model. 

PACS numbers; 74.25. Gz, 74.72. Jt, 78.30.Er 
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Since the last few years, identifying the symmetry of 
the pairing state has been the expected major step to- 
wards an understanding of high Tc superconductivity. 
The controversy between s-wave or d-wave pairing is not 
however resolved: some experiments appear to advocate 
strongly in favor of a d^^-y^ symmetry [Q , whereas others 
seem to show a significant s-wave contribution [Q . With 
respect to this problem, inelastic light scattering has been 
shown to be a powerful tool because besides probing the 
bulk (in contrast with photoemission and tunneling), the 
selection rules on the polarization of incoming and out- 
going light make the spectra sensitive to the wavevector 
of the electronic excitations (|-D . Usual theoretical ap- 
proaches of electronic Raman scattering ignore other ex- 
citations such as phonons. The experimental major diffi- 
culty lies however in the fact that in La2-xSrxCuOi, 
YBa2Cu307-5 (Y-123), Bi2Sr2CaCu20s+5 (Bi-2212) 
and Tl2Ba2CuOQ+s the phonons mask the electronic ex- 
citations, which hampers an accurate determination of 
their spectrum in the superconducting state 

In this letter, we report for the first time pure - no 
phonon structures superimposed - electronic Raman spec- 
tra.in HgBa2CaCu20e+s (Hg-1212) single crystals. This 
compound belongs to the new mercury family where the 
phonons do not mask the low energy electronic spectrum 
for electric fields within the planes [|l5|,^. Its crys- 
tallographic structure is purely tetragonal {^D^h) [0, 
which allows an unambiguous comparison with theoreti- 
cal calculations based on tetragonal symmetry. IIg-1212 
is therefore particularly well suited for the study of the 
superconducting gap. Our most striking result is that the 



low frequency behavior of the Big spectrum displays a 
quasi-linear dependence; this has not been reported sofar, 
presumably because the determination of the spectrum 
in this frequency range is usually hindered by phonons. 
Such an observation strongly suggests that the nodes are 
shifted from the [110] and [110] directions, a conclusion at 
odds with the generally accepted dx2_y2 gap symmetry. 
We believe that our new and unambiguous experimen- 
tal data should be taken into account in any theoretical 
approach of high- Tc superconductors. 

Raman spectra have been obtained from high-quality 
single crystals of IIg-1212. The crystals were grown by 
a single step synthesis as previously described for Hg- 
1223 |18|. They are perfect parallelepipeds with typical 
0.5x0.5 mm^ cross section and thickness 0.3 mm. The 
crystals were characterized by X-ray diffraction and en- 
ergy dispersive X-ray analysis [jl^. The [100] crystallo- 
graphic direction lies at 45° of the edge of the square 
and the [001] direction is normal to the surface. DC 
magnetic susceptibility measurements performed with a 
SQUID magnetometer yield = 126 ± IK. 

The measurements were performed with a standard 
Raman set-up using a single channel detection and the 
Ar'^ laser 514.52 nm line [|l^. The spectral resolution 
was set at 3 cm~^. The crystals were mounted in vac- 
uum (10^^ mbar) on the cold finger of a liquid helium 
cryostat. The temperature was controlled by a Si diode 
located inside the cold finger. The incident laser power 
density onto the crystal surface was kept below 15W/cm^ 
in order to avoid heating of the crystal during the low- 
temperature runs. The incidence angle of the incoming 
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light was 60° and the scattered Hght was collected along 
the normal to the crystal surface. The polarizations of 
the incident and scattered lights are denoted in the usual 
way. x: [100] (a axis), y: [010], z: [001] (c axis), x': 
[110], y': [110]. In order to compare our experimental 
data with theoretical calculations |5lP,|20|, the pure Aig 
(xx), i?2g (xy) and Big (x'y') symmetries are needed. 
This requires that the incident electric field lie within 
the xy plane. The incident angle being non zero, the 
crystal must be rotated by 45° after measuring in the 
Aig and B2g channels in order to get the Big spectrum. 
The impact of the laser beam onto the crystal surface was 
precisely located in order to probe the same crystal area 
for the three symmetries. Finally, a very weakly diffusive 
spot was carefully chosen for each crystal to minimize the 
amount of spurious elastic scattering. 

The raw Raman spectra at T=13 K of the Hg-1212 
single crystals for the Aig, Big and B2g symmetries are 
displayed in Fig. la). Remarkably and in sharp contrast 
with Y-123 and Bi-2212 l]9|,|l3|, no Raman active mode 
hinders the Big and B2g spectra. Weak peaks are seen 
at 388 and 576 cm~^ in the Aig spectrum but their en- 
ergy location is such that no correction is actually needed 
to discuss the electronic spectra [^ . Note that the resid- 
ual elastic scattering below 50 cm~^ is only seen in the 
Aig spectrum and is very weak. Therefore, we are in a 
position to turn immediately to the analysis of the data 
without dealing with the delicate handling of "phonon 
subtraction" . 

Two well marked maxima are observed around 530 
cm~^ and 800 cm~^ for the Aig and Big symmetry re- 
spectively. The ratio of the peak energies is ~ 1.5 while 
it was found -2 for Bi-2212 and Y-123 |Q. No clear 
maximum appears in the i?2g symmetry, for which actu- 
ally, the scattered intensity levels out smoothly around 
530 cm^^. The intrinsic (e.g. after subtracting the dark 
current of the photomultiplier) scattered intensities in 
Aig, Big and i?2g symmetries drop close to a zero count 
rate at zero frequency whereas, in Bi-2212, a trend to- 
wards finite intensity was claimed to be observed in Aig 
symmetry [^. We note that in Bi-2212, the presence of 
numerous phonon structures at low frequency especially 
in x'x' polarizations (xx in the notation of ref. makes 
it difficult to estimate quantitatively the electronic scat- 
tered intensity. In this respect, we believe that our data 
are more reliable. 

Fig. lb) displays the imaginary part Im[x(w, T=13K)] 
of the electronic response functions at T=13 K associ- 
ated to the Aig, Big and B2g symmetries. They were ob- 
tained from the raw spectra after subtracting the residual 
Rayleigh scattering, and correcting for the thermal Bose- 
Einstein factor 1 -I- n{u)) = [1 — exp{~huj / kBT)]~^ . The 
dashed lines in Fig. lb) are the imaginary part Im[x(w 
,T=150K)], of the response functions at T=150 K dis- 
played for clarity after smoothing the spectra. Note that 
the three response functions vanish at zero frequency. 



The main new observations which emerge from Fig. lb) 
are: i) the Big response function exhibits a clear decrease 
of the electronic scattering rate at low energy with re- 
spect to the normal state, unlike the Aig and the i?2g 
response functions where the difference is small, ii) the 
normal state intensity is recovered for all symmetries be- 
yond the maximum, iii) the two distinct maxima ob- 
served for both Aig and Big symmetries disappear above 
Tc- iiii) the scattered intensity in the Big channel ex- 
hibits a quasi-linear increase. 

Our experimental data are of special interest in the low 
frequency regime. Indeed the Big and B23 symmetries 
are the ones which are expected to have the most differ- 
ent behavior {u^ and uj respectively) in the d^^-yz model. 
Two fits have been performed for each symmetry: a first 
one to a power law and a second one to the simple 
polynomial function huj +cuj^. The a exponents calcu- 
lated below 300 cm^^ for Aig, B2g and Big symmetries 
are 1.3 ±0.3, 0.8 ±0.2, 1.5 ±0.5 respectively. The power 
law fit shows clearly that the low energy range of the Big 
spectrum does not display a uj^ dependence, in contrast 
to the claim of ref. []5|,p ,pOt . A linear u dependence in the 
low frequency regime is quite compatible with our three 
response functions. In the second fit the ratio T=b/ (cw^), 
which gives the relative weight of the linear and the cubic 
variation, is r = 2±1, 4±1, 3.3±0.7 and 3.4±0.5 at 300, 
400, 500 and 600 cm^^ respectively, much larger than an 
earlier estimate . This provides quantitative evidence 
for the predominance of the linear part compared to the 
cubic part especially for the Big symmetry. We suspect 
that the uj^ fit to previous experimental measurements is 
less reliable due to the phonon background PJT^. 

We now turn to the comparison of our data with ex- 
isting theories. In the limit where the superconducting 
coherence length is much smaller than the optical pene- 
tration depth, the imaginary part of the unscreened re- 
sponse function for T — > is given by [^ : 
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Np is the density of states for both spin orientations at 
the Fermi level, and the brackets indicate an average over 
the Fermi surface. stands for the superconducting, k- 
dependent gap. 7fe is the Raman vertex. 

Obviously, our results cannot be accounted for with an 
isotropic gap combined with a cylindrical Fermi surface. 
Several A^ and 7fc distributions have been proposed to 
reproduce the electronic excitations. Taking a Gaussian 
distribution of A^ 's |^,^, we obtain for the Aig and Big 
spectra the fits shown in Fig. lb). We have added a w^/^ 
function to the bracketed expression at energies greater 
than the A^ values, in order to mimic the asymptotic re- 
covery of the normal state behavior. The low energy part 
of the response function in Big symmetry is not satisfac- 
torily described by this type of calculation and moreover 
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this numerical fit does not yield any information on the 
gap symmetry. 

We have then attempted to fit our Aig, Big and B2g 
response functions with the calculations of Devereaux et 
al. for a rfa;2_y2 gap which were found to describe 
satisfactorily the Bi-2212 results. After adjusting the rel- 
ative energy location of the maxima in Aig and Big, we 
find a fair agreement with our data in the linear depen- 
dence of the low energy part of the Aig and i?2g sym- 
metries (see Fig. 2). However, the low frequency behav- 
ior of the theoretical Big spectrum is incompatible with 
our experimental results. This was already clear from 
the above discussion of the exponents and of the ratio of 
the linear and cubic contributions. We have also done 
the calculation for the simplest s-wave anisotropic gap 
A{9) = Aq + Al cos^ 26, not finding a satisfactory agree- 
ment because no hint of a threshold corresponding to the 
minimum gap is seen in our data. 

Let us therefore consider the implications of our exper- 
imental findings. We analyze first the Big and B2g sym- 
metries, which are easier to interpret due to simple sym- 
metry considerations and are not affected by screening 
in contrast with the Aig symmetry. In the Big case, the 
Raman vertex 7^ is zero by symmetry for = ky, hence 
the electronic scattering is insensitive to the gap structure 
around 45° ; in contrast, the ~ and ky k Q regions 
do contribute, giving weight to the gap Aq in these di- 
rections. Conversely, in the B^g case, the Raman vertex 
7fc is zero by symmetry for k^ = Q ov ky = Q, and non 
zero elsewhere, hence provides weight in the k^ — ky di- 
rection. In the dx'^-yi pairing state, these considerations 
imply that the Big symmetry is insensitive to the nodes 
at 45° (hence the uj'^ dependence) and displays a maxi- 
mum at 2Ao, whereas the B2g symmetry exhibits a linear 
frequency dependence (because it probes the nodes) and 
a smeared gap. An inescapable consequence of the linear 
low frequency dependence is the existence of a density 
of states which increases linearly with frequency, most 
naturally due to nodes in the gap. Since these nodes are 
probed in both Big and B2g symmetries, they cannot he 
(only) located in the 45° direction. 

It could be argued that impurities are responsible for 
the observed low frequency density of states in the Big 
channel. Devereaux p2| has shown that, for a d^2_y2 
gap, impurities induce for Big a linear rise of the elec- 
tronic scattering at low frequency, crossing over to the 
dependence at higher frequency. However the lin- 
ear behavior in Big extends in our experiment up to 600 
cm"-'^, yielding a cross-over energy lj* (FA)^/^ of order 
600 cm~^. In the most favorable case of unitary scatter- 
ers ||2^ , this would imply a scattering rate F of order 
of the gap A itself. Such a high scattering rate should 
strongly reduce the critical temperature, in contradic- 
tion with the optimal 126 K critical temperature of our 
samples. Therefore this explanation looks unlikely. We 
believe that our data are representative of a pure com- 



pound and that the linear contribution found in the Big 
spectrum is intrinsic. 

We are thus left with the conclusion that the nodes are 
shifted from the 45° direction. Hg-1212 being tetragonal, 
we cannot ascribe this to an orthorhombic distorsion as 
could be the case in Y-123. We therefore explore this shift 
more quantitatively. As a first approach, we take the sim- 
ple model 7_Big ^ cos 20 , 7523 ~ sin 26* for the Raman 
vertices [^,^ and a cylindrical Fermi surface. We start by 
using a "toy model" where we compute the Big and i?2g 
spectra for an order parameter A{9) = Aq cos(20 — 2a) 
obtained by artificially rotating the d^^^-y^ one by an an- 
gle a.. This is intended to check the sensitivity of the 
Raman spectra to the nodes location (note that a 22.5° 
rotation would naturally make Big and B2g identical). 
We find that a rotation by a «10° is barely noticeable in 
the spectra : the linear rise at low u produced in Big by 
this rotation is comparable to the experimental accuracy. 
This suggests that the shift of the nodes away from 45° 
is quite sizeable. 

Because of the tetragonal symmetry, one node located 
at 6^0 ( 7^ 45° or 0° ) implies 8 nodes at ±6*0 -f n7r/2 (n 
= 0,1,2,3 ). A very simple model corresponding to this 
situation is A{9) = Aq [cos (46') + s] with < s < 1, an 
order parameter which has the Aig symmetry p^ ]. The 
gap has a maximum Ao(H- s) for 9 = Q, seen in Big, and 
a secondary maximum Aq (1 — s) for 9 = 7r/4, seen in S23, 
while the node lies at 6*0 = (1/4) arccos(-s). This model 
could also be loosely called a g-wave model or a (super) 
extended s-wave model (in contrast to the d-wave model, 
which has the Big symmetry). With one more parame- 
ter we might decouple the position of the nodes and the 
size of the gap maxima. We choose the parameter s to 
account for the relative peak positions in Big and i?2g ( 
inasmuch as we consider that i?2g has a very broad peak 
around 500 cm~^ ). It can be seen in Fig. 2 that the Big 
and B2g spectra calculated in this framework for s=0.2 
(giving 00 = 25°) are in good agreement with experi- 
ment for the low frequency behavior. The discrepancy 
at high frequency, most conspicuous in B2g, is assigned 
to finite electronic scattering rate in the normal state. 
This should be also accounted for, but this is beyond the 
scope of this letter. In order to remove the singularities 
at the gap maxima, a smearing function with width pro- 
portional to frequency [2^ has been incorporated in the 
calculation. 

Our model explains the different positions of the Big 
and B2g peaks (also found in Bi-2212 ) by the gap 
anisotropy. This leads us to try and connect the posi- 
tion of the peak found in Aig symmetry to the other 
peaks. The Aig symmetry is screened and the peak posi- 
tion can not be related directly to A{9). Since a constant 
vertex is completely screened, the simplest model giving 
a nonzero result for Aig is 7Aig ^ cos 40. However, in 
our one parameter model, the position of the nodes is 
imposed by the relative peak position in Big and B2g ■ 
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This leads to the 'yAig cos 40 being essentially zero at 
the nodes, and yields a very small low u linear contribu- 
tion in contrast with experiment. We have thus included 
higher harmonics. The calculation shown in Fig. 2 is done 
with jAig ~ cos 46 — 3 cos 89 . The agreement is quite 
satisfactory both with respect to the low to behaviour and 
the peak position. However it is clear that this could even 
be improved in a more complicated two parameter model 
for A{9), decoupling the node position and the gap max- 
ima, and also including a dependence for the density of 
states. 

In conclusion, we have presented for the first time, pure 
electronic Raman spectra of Hg-1212 single crystals. Our 
most significant result is the observation of an intrinsic 
linear co dependence -not only in the i?2g spectrum- but 
also in the Big spectrum. This is an important clue in 
order to locate the nodes, which was not reported sofar 
in other compounds presumably because the experimen- 
tal determination of the low frequency regime is usually 
hampered by phonons. Our observations advocate in fa- 
vor of a shift of the nodes from the [110] and [110] direc- 
tions, which is inconsistent with the simple d^2_y2 model. 

We thank M. Cyrot and E.Ya. Sherman for very fruit- 
ful discussions. 
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FIG. 1. a) Raw Raman spectra of the Hg-1212 single crys- 
tals in Aig, Big and B2g symmetries at T=13 K b) Imaginary 
parts of the response function lm[x{LL>,T = 13iir)] in the Aig, 
Big and B2g. The dashed lines represent Im[x(t<-', T = 150K)]. 
The solid line is a fit to the Aig and Big spectra of a Gaus- 
sian distribution of gaps with 2Ao — 360 cm~^ and a — 200 
cm^^ for Aig and 2Ao=670 cm^^ and a =180 cm^^ for Big 
(Ao and a are the mean and the standard deviation of the 
Gaussian respectively) 

FIG. 2. From left to right: Aig, Big and B2g data and cal- 
culations for the dx2-y2 model (dashed line) and for a model 
with 8 nodes lying at 9o = 25° (full line) 
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